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AEROSOL  CHARGE  DISTRIBUTIONS  PRODUCED 
BY  RADIOACTIVE  IONIZERS 


1.  INTRODUCTION 

The  mechanical  generation  of  particles  often  produces  an  aerosol  that  is  highly  charged.  The 
charge  can  be  neutralized  by  passing  the  aerosol  through  an  irradiated  region  where  the  production  of 
ions  creates  a  zone  of  high  electrical  conductivity.  In  this  region  the  aerosol  charge  is  conducted  from 
the  particles  to  the  walls  of  the  neutralizer.  The  charge  distribution  on  an  aerosol  that  has  undergone  a 
neutralization  process  is  often  assumed  to  be  described  by  Boltzmann  statistics  with  no  net  charge 
(equal  numbers  of  positively  and  negatively  charged  particles).  While  this  description  of  charge  neu¬ 
tralization  is  adequate  for  many  processes  in  which  the  terminal  charge  distribution  is  of  little  conse¬ 
quence,  it  is  not  adequate  for  certain  applications.  The  differential  mobility  analyzer  (DMA)  size 
spectrometer,  in  particular,  requires  accurate  knowledge  of  the  charge  distribution  in  order  to  invert 
the  mobility  distribution  measured  with  the  DMA  to  a  size  distribution. 

Undoubtedly  the  most  widely  used  theory  to  predict  the  charge  distribution  is  Boltzmann  statis¬ 
tics  (Keefe,  Nolan,  and  Rich,  1959).  This  theory  does  not  allow  for  the  differences  in  positive  and 
negative  ions  and  thus  predicts  symmetric  charging;  it  also  does  not  give  the  time  required  to  achieve 
equilibrium.  Diffusion  charging  theories  allow  for  differences  in  the  diffusion  and  mobility  coeffi¬ 
cients  of  the  positive  and  negative  ions  and  predict  that  charging  in  a  balanced  bipolar  ion  environ¬ 
ment  will  produce  a  nonsymmetric  aerosol  charge.  (For  a  review  of  bipolar  charging  theories,  see 
Hoppel  and  Frick,  1986.)  Since  the  charging  depends  on  the  ionic  environment,  an  important  first 
step  is  to  understand  the  processes  that  occur  in  the  ionizer.  This  study  was  motivated  by  confusing 
experimental  results  obtained  when  we  attempted  to  test  the  validity  of  different  charging  theories  by 
measuring  the  ratio  of  positively  to  negatively  charged  particles  exiting  neutralizers.  It  soon  became 
apparent  that  the  measured  charge  ratio  was  as  much  a  function  of  the  ionizer  design  as  it  was  of  the 
charging  theory.  This  led  us  to  study  the  factors  that  determine  the  ion  concentrations  in  an  ionizer. 
Much  of  our  original  confusion  stemmed  from  our  ignorance  of  what  was  actually  occurring  in  and 
downstream  of  the  ionized  region. 

This  report  is  a  study  of  the  physical  mechanisms  that  control  the  aerosol  charge  distribution 
established  by  radioactive  neutralizers.  The  computations  are  of  two  types:  (a)  modeling  the  ion 
environment  as  the  air  moves  over  and  downwind  of  a  region  of  high  ionization,  and  (b)  calculating 
the  aerosol  charge  distribution  in  ion  environments  with  various  polar  ion  ratios.  Measurements  were 
made  of  (a)  the  polar  ion  concentrations  downwind  of  radioactive  sources  and  (b)  the  degree  of  asym¬ 
metric  charging  produced  by  several  ionizers,  including  a  prototype  neutralizer  that  produces  a  nearly 
symmetrical  equilibrium  charge  distribution  by  adjusting  the  ratio  of  the  polar  ion  concentrations  to 
compensate  for  the  differences  in  the  ion  properties. 
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2.  BIPOLAR  CHARGING  IN  AN  UNBOUNDED  REGION 

The  statistical  charge  state  of  a  monodisperse  aerosol  is  determined  by  the  following  system  of 
balance  equations  (see,  for  example,  Hoppel  and  Frick,  1986). 

=  q  -  cmxn2  -  n/j:  -  nt  £  (i$Np ,  (1) 

at  k  =o  *  =»  1 

=  q  -  anin2  -  n2  flgWJ0  -  n2  £  (2) 

it  =o  *  =  i 

~  =  nfilW  +  n&M"  -  "itio'No  -  »20$ty>  •  (3) 


=  n1/31V)N1<*"1)  -  n  1/31^W1<*)  +  n20#  +  lW,(*  +  l)  -  n2^N\k)  - 


=  n2/3$-1)Nf _1>  -  n2(3&)N?)  +  Mt2  +  1)^  +  1)  -  , 


Z  =  N0  +  £  N,(t>  +  £  Af  >  ,  (6) 

Jt  =  i  t  =  i 

where  n,  and  n2  are  the  positive  and  negative  ion  concentrations  respectively,  j3,y  )  are  the  attachment 
coefficients  between  an  ion  of  polarity  i  and  a  particle  with  k  charges  of  polarity  j ,  Nfk)  is  the  con¬ 
centration  of  positive  or  negative  particles  with  k  charges,  Z  is  the  total  concentration  that  is  held 
constant  with  time,  q  is  the  ionization  rate,  and  a  is  the  recombination  coefficient.  There  is  one 
more  equation  than  the  number  of  variables.  In  the  limit  as  the  system  approaches  steady  state,  Eq. 
(3)  is  not  independent  of  Eqs.  (4)  and  (5);  we  therefore  chose  to  eliminate  Eq.  (3).  For  symmetric 
charging,  the  concentrations  of  positive  and  negative  ions  are  the  same,  as  are  the  numbers  of  posi¬ 
tively  and  negatively  charged  particles;  thus  charge  conservation  is  inherent.  For  the  nonsymmetric 
case,  these  equations  need  to  be  supplemented  with  the  charge  conservation  equation: 

P  p  nn 

n,  -  n2  +  £kN\k)  -  =  —  ,  (7) 

*=i  k= l  e 

where  p0  is  the  initial  space  charge  (obtained  from  the  initial  conditions). 

It  is  often  convenient  to  simplify  the  equations  by  assuming  fixed  ion  ratios  and  steady  state  con¬ 
ditions,  for  which  case  the  aerosol  balance  equations  decouple  from  the  ion  balance  equations.  The 
principle  of  detailed  balancing  (Gunn,  1955)  requires  that 

=  njffiNp.  (8) 

In  Eqs.  (4)  and  (5)  this  requires  that  the  first  and  fourth  terms  and  the  second  and  third  terms  balance 
separately.  (The  necessity  of  this  requirement  at  steady  state  can  be  demonstrated  by  terminating  the 
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above  system  of  equations  at  some  finite  value  of  k.)  Solving  for  the  number  of  particles  with  k 
charges  in  terms  of  uncharged  particles  gives 


N,(k) 


No 


n  o(m) 

^  II 


m  =  1 


(9) 


The  charge  distribution  //*  *  describes  the  fraction  of  particles  carrying  k  charges  of  sign  i . 
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where  the  numerator  is  Eq.  (9)  and  the  denominator  is  obtained  from  Eqs.  (6)  and  (9): 
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These  equations  do  not  include  any  diffusive  or  electrically  driven  fluxes  of  ions  or  aerosols  to 
the  walls  of  the  container.  Because  of  the  higher  diffusion  coefficient  of  ions,  their  diffusive  loss  to 
the  walls  is  much  greater  than  for  the  aerosol.  The  loss  of  ions  to  the  walls  is  considered  in  Section 
3. 

The  most  widely  used  diffusion  theory  for  determining  the  ion-aerosol  attachment  coefficients  is 
that  of  Fuchs  (1964).  Fuchs’  theory  is  reviewed,  and  a  theory  that  includes  three-body  trapping  is 
presented  in  a  paper  by  Hoppel  and  Frick  (1986),  which  provides  attachment  coefficients  for  ions  of 
mobilities  of  1.2  and  1.35  cm2  V-1  s-1  (atmospheric  ions  at  standard  temperature  and  pressure 
(STP),  (Mohnen,  1977)).  The  use  of  these  values  of  mobility  together  with  the  Nolan  value  of  1.4 
x  10~6  cm3  s'1  were  motivated  by  our  interest  in  the  interaction  of  ions  and  aerosols  in  the  atmo¬ 
sphere  (Hoppel,  1985).  For  higher  ion  concentrations  (shorter  ion  lifetimes)  as  exist  downwind  of  the 
ionizer,  more  appropriate  values  of  positive  and  negative  ion  mobilities  would  be  1.3  and  1.8 
cm2  V-t  s'1  respectively,  at  STP  and  an  ion-ion  recombination  coefficient  of  1.6  x  10-6  cm3  s_1. 
The  resulting  negative  to  positive  mobility  ratio  of  about  1 .4  is  in  good  agreement  with  the  measure¬ 
ments  downwind  of  the  ionizers  reported  in  Appendix  B.  The  recombination  coefficients  have  been 
recalculated  with  these  values  of  mobilities  and  are  given  in  Appendix  D,  together  with  a  discussion 
of  the  method  of  calculations.  In  this  report  we  use  the  values  of  j3  given  in  Appendix  D. 

2.1  Symmetric  Charging 


Symmetric  charging  is  an  ideal  process  for  which  the  positive  and  negative  ions  are  assumed  to 
be  identical  except  for  their  polarity,  thus  $tJ  =  (5j, ,  /3„  =  0Jjy  and  n,  =  n2-  Consequently,  the 
charge  distribution  is  symmetric  and  Eqs.  (9)  and  (10)  become 

k 


z 


m  - 1 


1  +  2 


p 

z 


n  0% 


(m  -  1 ) 


(12) 


k  =  I  m  -  I 


Figure  1  shows  the  numerical  solution  (discussed  later)  to  Eqs.  (5)  through  (II)  for  the  case  when 
q  —  104  cm"3  s-1.  The  calculations  are  for  0.05-^im  radius  particles,  using  values  of  0  for  ions  of 
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Fig.  1  —  Time-dependent  symmetric  charging  of  an  initially  uncharged  monodisperse 
aerosol  (r  =  0.05  nm  )  by  ions  of  mobility  /t  =  1.3  cm :  V”1  s' 1 

mobility  1.3  cm2  V-1  s_l.  If  the  calculations  are  repeated,  using  values  of  0  for  1.8  cm2  V-1  s~‘ 
for  both  positive  and  negative  ions,  the  equilibrium  charge  distribution  is  nearly  identical  but  is 
achieved  somewhat  sooner.  The  equilibrium  distribution  is  shown  in  Fig.  2  as  a  function  of  size.  The 
equilibrium  charge  distribution  is  nearly  independent  of  the  ionic  properties  in  the  symmetric  case.  It 
is  interesting  to  note  that  the  Boltzmann  theory  predicts  that  the  charge  distribution  (shown  by  the 
dashed  line  in  Fig.  2)  is  totally  independent  of  the  ionic  properties.  The  fact  that  the  charge  distribu¬ 
tion  is  nearly  independent  of  ionic  properties  for  symmetric  charging  will  later  be  used  as  an  argu¬ 
ment  for  the  design  of  a  neutralizer  that  produces  a  near-symmetric  aerosol  charge  distribution. 

2.2  High  Ionization  Rate— Asymmetric  Charging 

In  a  region  of  very  high  ionization  (and  aerosol  concentrations  limited  so  that  n  »  Z),  the  ion 
densities  are  determined  by  the  balance  between  the  rates  of  ionization  and  recombination.  These  are 
symmetric  processes  and  produce  an  environment  where  the  polar  ion  concentrations  are  equal.  In  the 
high-ionization  limit,  any  ionic  imbalance  caused  by  the  unequal  attachment  of  positive  and  negative 
ions  to  the  aerosol  has  negligible  effect  on  the  ion  concentrations.  In  this  case  Eqs.  (9),  (10),  and 
(11)  together  with  the  assumption  of  equal  ion  densities  can  be  used  to  calculate  the  asymmetric 
charge  distribution.  Assuming  values  of  0  corresponding  to  positive  ions  of  mobility  1.3  and  values  of 
0  for  negative  ions  of  mobility  of  1.8  (given  in  Appendix  D),  we  obtain  Fig.  3.  This  charge  distribu¬ 
tion  is  obviously  much  different  from  that  shown  in  Fig.  2  and  illustrates  that,  for  real  situations 
where  aerosols  are  neutralized  in  a  region  of  high  ionization,  we  should  expect  the  charge  distribution 
to  be  asymmetric. 
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Figure  4  displays  the  charge  distribution  for  particles  of  radius  0.01,  0.05,  and  0.2  /im  as  a 
function  of  the  number  of  elementary  charges.  The  dashed  line  gives  the  charge  distribution  derived 
by  Gunn  (1955)  for  the  case  when  the  radius  is  larger  than  0.05  /im.  The  Gunn  expression  is 

N(k)  e  f  ~[k  ~  {rKT/e2)ln(jilnl/n2n2)]2 

~Z~  ~  (2irrKT)1/2  l  2^T 

where  K  is  Boltzmann’s  constant,  T  is  the  absolute  temperature,  and  m  and  n2  are  the  ionic  mobili¬ 
ties.  Gunn’s  expression  was  developed  for  particles  of  radii  greater  than  0.05  /tm;  this  creates  an 
inconsistency  regarding  the  sum  of  charged  and  uncharged  particles  for  smaller  sizes,  as  is  evident  for 
the  0.01 -/tm  size.  Gunn’s  expression  is  similar  to  the  Boltzmann  charge  distribution  except  that  it  is 
offset  from  k  =  0  by  a  factor  related  to  the  ratio  of  the  polar  conductivities.  For  the  case  shown 
(high  ionization)  we  have  assumed  that  the  ion  concentrations  are  equal,  and  therefore  the  offset  is 
due  only  to  the  difference  in  the  polar  ion  mobilities. 


CHARGE 

Fig.  4  —  Fraction  of  particles  carrying  k  elementary  charges  for  particles  of  0.01,  0.05, 
and  0.1  radius  calculated  for  positive  and  negative  ion  mobilities  of  1.3  and  1.8 
cm2  V"1  s-1  respectively.  The  dashed  line  is  from  Eq.  (13). 


Once  the  sample  leaves  the  region  of  high  ionization,  the  ion  concentration  decays  very  rapidly, 
and  it  may  be  impossible  to  maintain  conditions  such  that  the  ion  concentration  is  determined  only  by 
symmetric  processes  (e.g.,  recombination).  An  imbalance  in  the  polar  ion  concentrations  develops 
when  the  ion  loss  to  aerosols  and/or  the  walls  is  comparable  to  loss  by  ion-ion  recombination. 

2.3  Low  Ionization  Rate  or  High  Aerosol  Concentrations,  Z  »  n 

For  a  low  ionization  rate  and/or  high  aerosol  concentration,  the  loss  of  ions  to  the  aerosol  parti¬ 
cles  is  the  dominant  sink  for  ions  (the  particles  act  as  recombination  centers  for  the  ions).  If  there  is 
initially  no  net  space  charge,  then  charge  conservation  requires  that  the  concentrations  of  positively 
and  negatively  charged  particles  be  equal.  Production  of  this  symmetric  aerosc!  charge  requires  the 


NRL  REPORT  9108 


positive  to  negative  ion  ratio  to  be  greater  than  unity  (Eq.  (9)).  This  imbalance  of  charge  on  the  less 
numerous  ions  must,  of  course,  be  balanced  by  a  small  difference  in  the  polar  charge  ratio  carried  by 
particles  (the  ratio  will  still  be  near  unity).  An  important  example  of  the  low  ionization  case  is  a  pol¬ 
luted  atmosphere  for  which  the  concentration  of  ions  due  to  natural  radioactivity  is  much  less  than  the 
concentration  of  aerosols.  In  a  neutralizer,  charge  equilibrium  must  be  achieved  rapidly,  which 
requires  the  high  ion  concentrations  provided  by  high  ionization  rates.  The  convergence  of  the 
charged  particle  ratio  to  unity  as  the  ionization  is  decreased  is  shown  in  Hoppel  and  Frick  (1986). 
This  calculation  has  been  replotted  in  Fig.  5  (using  the  values  of  <3  given  in  Appendix  D)  in  terms  of 
the  ratio  of  ion  loss  by  recombination  to  the  loss  by  attachment  of  aerosols. 


\  qa 


Fig.  5  —  Imbalance  in  positive  and  negative  particles  as  a  function  of  the  dimension¬ 
less  parameter  (qa)'n/0Z  for  particles  of  0.01  and  0.1  jtm  radius 


2.4  Relaxation  Time  for  Aerosol  Charging 

The  time-dependent  charging  and/or  discharging  of  aerosols  has  been  modeled  by  applying  a 
Runga-Kutta  method  to  solve  Eqs.  (1)  through  (7),  for  which  the  values  of  (3  were  taken  from 
Appendix  D.  Solutions  for  particles  of  radius  0.01  and  0.1  n m  are  shown  in  Figs.  6  and  7  and  illus¬ 
trate  the  approach  to  equilibrium  as  an  initially  uncharged  aerosol  is  exposed  to  a  region  where  the 
ionization  rate  is  105  ion  pair  cm-3  and  the  recombination  coefficient  is  taken  to  be  1.6 
x  10~6  cm3  s-1.  The  numbers  on  the  curves  correspond  to  the  number  of  elementary  charges  on  the 
particles.  The  total  particle  concentration  in  each  case  is  104.  For  0.01-^m  particles,  there  are  no 
multiply  charged  particles.  If  we  were  to  assume  that  the  ionization  rate  is  sufficiently  high  that  ion 
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recombination  dominates  ion  loss  to  the  particles  and  the  charging  is  symmetric,  then  the  ion  equa¬ 
tions  decouple  from  the  aerosol  equations  and  an  analytic  solution  is  possible.  The  solutions  for  the 
ion  concentrations  and  fraction  of  singly  charged  particles  are: 


n  +  =  «  _  =  (qa)U2  tanh  [(a^)1/2t] 


0o 

20o  +  0, 


cosh  (2 t'faq  )  +  1 


20o  +  0, 


These  analytic  solutions  are  shown  by  dashed  lines  in  Figs.  6  and  7  and  give  nearly  the  same 
result  as  the  rigorous  solution  for  the  0.01-gm  size.  For  the  0.1-gm  particles  (Fig.  7),  multiply 
charged  particles  are  important.  Since  a  k  +  1  charged  particle  is  derived  from  a  particle  with  k 
charges,  it  appears  at  a  later  time.  The  dashed  lines  are  again  obtained  from  Eqs.  (14)  and  (15)  where 
the  values  of  0  used  were  the  coefficients  for  single  charging  and  discharging  events;  multiply 
charged  particles  are  neglected.  Equations  (14)  and  (15)  are  given  here  so  that  estimates  of  the  time 
required  for  equilibrium  over  a  radioactive  source  can  be  obtained  without  resorting  to  the  numerical 
solution  of  the  complete  set  of  equations. 

Appendix  C  presents  further  graphical  solutions  of  aerosol  charging  and  discharging  under  vari¬ 
ous  conditions. 

2.5  Aerosol  Charging  After  the  Ionization  Is  Terminated 

In  the  event  that  the  ionization  source  is  abruptly  removed  after  charge  equilibrium  is  achieved, 
then  the  particles  will  age  in  a  region  of  rapidly  decaying  ion  concentrations.  When  the  positive  and 
negative  ions  have  identical  properties,  the  charge  on  the  aerosols  is  symmetric  and  remains  sym¬ 
metric  in  the  region  of  decaying  ion  concentrations  because  the  polar  ion  ratio  remains  at  unity.  How¬ 
ever,  if  for  some  reason  there  is  a  slight  imbalance  in  the  number  of  positive  and  negative  ions,  then 
conservation  of  charge  requires  that  there  be  an  ever-increasing  ion  ratio  as  the  ions  decay.  In  this 
case  the  aerosol  charge  distribution  changes  as  a  result  of  the  changing  ion  ratio,  and  equilibrium  will 
not  exist.  The  initial  ion  imbalance  may  be  caused  by  more  rapid  diffusion  of  one  ion  species  to  the 
walls  or  by  some  more  exotic  mechanism  such  as  the  expulsion  of  secondary  electrons  during  radioac¬ 
tive  decay.  Figure  8  illustrates  how  the  changing  ion  ratio  during  decay  would  affect  the  aerosol 
charge  distribution  if  the  initial  ion  ratio  were  0.98.  In  the  realistic  case  of  ions  with  unequal  diffusion 
coefficients,  there  will  always  be  unequal  wall  losses  that  will  be  magnified  in  the  ionic  ratio 
downwind  of  the  ionizing  source,  where  the  ion  concentrations  are  decaying  because  of  recombina¬ 
tion.  This  indicates  that  the  equilibrium  charge  established  in  the  ionizer  will  be  modified  after  the 
aerosol  has  been  removed  from  the  source  region.  Since  in  most  ionizers  there  will  be  a  region  of 
decaying  ion  concentrations  downwind  of  the  radioactive  source,  either  in  the  ionizer  itself  or  in  the 
plumbing,  the  ability  of  an  ionizer  to  produce  an  equilibrium  charge  distribution  on  an  aerosol  is 
brought  to  question. 

3.  ION  CONCENTRATIONS  AND  AEROSOL  CHARGING  DURING  FLOW  THROUGH  A 
CHANNEL  CONTAINING  AN  IONIZER 

There  are  always  ion  losses  to  the  walls  and  the  plumbing  between  the  ionization  source  and  the 
measuring  device.  To  explore  the  possibility  of  modification  of  the  aerosol  charge  distribution 
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Fig.  8  —  The  increase  in  the  ratio  of  the  polar  ion  densities  caused  by  an  initial  ion 
imbalance  (2%)  and  its  effect  on  the  charge  distribution  of  0.01  gm  radius  particles 


downwind  of  the  ionization  region,  we  have  modeled  the  ion  concentrations  in  a  channel.  The  equa¬ 
tions  that  govern  the  ion  concentrations  are: 


V  7,  =  q  -  an  |/i2 . 


(16) 


and 


V  ■  J 2  =  q  -  an  {r,2 , 


V  •  E  =  —(n  |  ~  n-,) , 

t 

and  the  ion  fluxes  are 


J  |2  —  i  En  1  2  ^  1,2  D  |  2^t  i  2  * 


(17) 


(18) 


(19) 


where  D\  2  are  the  diffusion  coefficients,  E  is  the  electric  field  resulting  from  any  ionic  space  charge, 
e  is  the  elementary  charge,  e  is  the  permittivity  of  air,  and  v  is  the  flow  velocity. 

These  equations  include  the  following  ionic  effects:  diffusion,  electrical  conduction  to  the  walls, 
formation  by  volume  ionization,  and  recombination.  Diffusion  and  conduction  parallel  to  the  flow  are 
assumed  to  be  negligible  compared  to  the  perpendicular  components,  and  the  flow  is  assumed  to  have 
a  parabolic  profile.  Also,  ion-aerosol  interactions  are  assumed  to  be  unimportant  in  establishing  the 
ion  concentrations;  an  assumption  that  is  only  valid  in  a  region  where  the  ion  concentrations  are  much 
larger  than  the  aerosol  concentrations. 
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3.1  Relative  Importance  of  Different  Processes 

For  illustrative  purposes  we  have  modeled  the  ion  concentrations  in  a  2-cm-wide  channel.  The 
air  flows  over  a  radioactive  source  where  it  reaches  an  equilibrium  ion  concentration  of  1  x  105. 
Figure  9(a)  shows  the  ion  profiles  at  various  distances  downstream  of  the  source  for  a  parabolic  flow 
profile  with  mean  velocity  of  6  cm/s.  The  solid  line  is  for  positive  ions,  and  the  dashed  line  is  for 
negative  ions  (of  mobilities  1.3  and  1.8  cm2  V-1  s_1),  respectively.  To  understand  the  relative  impor¬ 
tance  of  the  various  physical  processes,  we  compare  modeled  results  for  which  a  specific  process  has 
been  disabled.  Figure  9(a)  shows  the  complete  solution  of  Eqs.  (16)  through  (19)  for  the  ion  profiles. 
The  difference  in  the  positive  and  negative  profiles  is  due  to  the  more  rapid  diffusion  of  negative  ions 
to  the  walls. 


ION  DENSITY  (x  105  cm"3) 


Fig.  9(a)  —  Complete  solution  for  the  evolution  of  the  ion  profiles  in  a 
2-cm  channel  downwind  of  a  bipolar  ion  source.  Solid  (dashed)  lines 
are  for  positive  (negative)  ions. 


Figure  9(b)  shows  that  if  ion  recombination  is  omitted  (a  =  0),  the  ion  concentration  decay  is 
greatly  slowed  as  the  air  flows  down  the  tube.  Since  recombination  is  proportional  to  the  product  of 
the  ion  concentrations,  recombination  is  the  dominant  loss  mechanism  at  high  ion  concentrations.  The 
profile  shown  in  Figure  9(b)  is  largely  the  result  of  ionic  diffusion  to  the  walls.  However,  the  electric 
field  caused  by  the  ionic  space  charge  is  also  important  in  that  it  attracts  negative  ions  from  the  walls 
and  repels  positive  ions  outward  toward  the  walls.  This  limits  the  degree  of  ionic  imbalance  that  can 
develop  because  of  diffusion. 

I  If  diffusion  to  the  walls  is  turned  off,  as  illustrated  by  Fig.  9(c),  the  concentrations  of  positive 

and  negative  ions  remain  equal  and  there  is  no  electric  field  to  influence  ion  migration.  The  concen¬ 
tration  profile  develops  because  of  the  parabolic  velocity  profile;  the  air  near  the  walls  flows  more 
slowly  and  the  ions  have  more  time  to  recombine.  (If  plug  flow  is  assumed,  then  the  concentration 
profile  remains  constant  with  height.) 
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Figure  9(d)  shows  that  if  electrical  conduction  is  turned  off  by  setting  the  ionic  mobilities  to 
zero,  diffusion  of  negative  ions  to  the  walls  exceeds  that  of  the  positive  ions,  producing  a  larger  polar 
ion  ratio  ( n  +  /n_ )  than  would  otherwise  occur.  This  effect  may  appear  to  be  a  small  contribution  to 
the  overall  ion  concentration,  but  its  effect  on  the  ion  ratio  is  quite  significant.  In  the  case  shown  the 
predominance  of  positive  ions  has  increased  from  8%  to  17%  (profile  average).  As  shown  earlier,  the 
ratio  of  the  polar  ion  densities  is  important  in  determining  the  ratios  of  positively  to  negatively 
charged  particles.  Since  the  electric  field  is  the  integral  of  the  space  charge,  the  effect  of  the  electric 
field  is  strongly  dependent  on  the  size  of  the  neutralizer. 


ION  DENSITY  (x  10s  cm'J) 


Fig.  9(d)  —  Change  in  solution  (a)  obtained  by  elimination  of  the 
mobilities  of  the  ions  in  the  electric  field  (n+._  =  0) 


3.2  Comparison  Between  Experiments  and  Modeling  Results 

Polar  ion  concentrations  were  measured  for  several  distances  downwind  of  a  radioactive  source 
(5.3  gCi  Am-241)  within  a  rectangular  (1.9  x  4.5  cm)  tube.  Ion  concentrations  were  calculated  from 
the  expression: 


where  i  is  the  ion  current,  e  is  the  electron  charge,  and  is  the  volumetric  flow  rate.  This  equation 
assumes  that  the  voltage  on  the  collector  is  large  enough  to  collect  all  the  ions,  and  therefore  the 
measured  concentrations  represent  the  mean  ion  concentrations. 

Figure  10  illustrates  an  experiment  in  which  the  radioactive  source  was  positioned  distances  of 
6,  12,  and  18  cm  from  the  shield  wires.  Shorter  distances  were  not  used  because  a  series  of  I-V 
curves  taken  during  a  preliminary  portion  of  this  experiment  revealed  that  current  saturation  was  diffi¬ 
cult  to  achieve  at  less  than  6  cm  (see  Appendix  B).  Ion  concentrations  at  the  6-cm  location  were  used 
as  an  initial  condition  in  the  model,  and  an  initial  concentration  profile  was  adopted  to  simulate  the 
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Fig.  10  —  Experimental  arrangement  for  measurement  of 
ion  concentrations  downwind  of  a  radioactive  source 


first  6  cm  of  travel.  Figure  11  shows  the  computed  concentrations  of  positive  and  negative  ions  12 
and  18  cm  downwind  of  the  source  (the  upper  set  of  dashed  lines),  compared  with  the  experimental 
data  (the  upper  solid  lines).  Closer  agreement  could  be  obtained  by  adjusting  the  values  used  for 
recombination  and  mobility;  however,  since  the  model  assumes  parallel  plates,  whereas  the  experi¬ 
mental  setup  has  side  walls,  precise  agreement  was  not  expected.  A  second  comparison  was  per¬ 
formed  in  which  the  6  cm  immediately  in  front  of  the  shield  wires  was  spanned  by  a  set  of  five  paral¬ 
lel  plates  (each  4.0  x  6.0  x  0.05  cm),  spaced  0.3  cm  apart;  the  source  was  located  6  cm  upwind  of 
this  “diffusion  cell”  (12  cm  from  the  shield  wires).  For  this  case,  the  measured  concentration  of  the 
positive  and  negative  ions  is  indicated  by  the  lower  solid  lines  of  Fig.  11,  and  the  computed  concen¬ 
trations  are  indicated  by  the  lower  set  of  dashed  lines.  In  addition  to  the  large  diffusive  losses 
expected  in  the  diffusion  cell,  there  is  an  increase  in  the  polar  ion  ratio.  This  is  because  the  smaller 
channels  of  the  diffusion  cell  reduce  the  electric  field  caused  by  the  ionic  space  charge,  therefore  a 
larger  polar  ion  ratio  can  develop. 

4.  SUMMARY  STATEMENTS 

The  following  statements  based  on  the  preceding  two  sections  summarize  the  mechanisms  that 
are  important  in  understanding  the  operation  and  design  of  radioactive  neutralizers. 

1.  In  a  region  of  high  volume  ionization,  the  polar  ion  concentrations  will  be  nearly  equal. 

2.  When  the  ion  concentrations  are  equal,  there  will  be  pronounced  asymmetric  charging  owing  to 
the  difference  in  the  positive  and  negative  ion  diffusivities  (mobilities). 

3.  In  the  region  of  decaying  ion  concentrations  downwind  of  a  source,  the  polar  ion  ratio  will  be 
changing  and  no  true  equilibrium  charge  distribution  can  be  maintained.  In  an  ionizer  this  is  due 
primarily  to  the  unequal  diffusion  of  positive  and  negative  ions  to  the  walls.  However,  even  in 
the  absence  of  walls,  the  charge  distribution  will  change  as  the  system  evolves  from  an  ion- 
dominated  regime  to  an  aerosol-dominated  one,  as  is  shown  in  Fig.  5. 

4.  Since  the  diffusion  coefficients  of  positive  and  negative  ions  are  unequal,  large  ratios  of  ion  den¬ 
sities  can  be  obtained  in  a  diffusion-dominated  region  downwind  of  the  region  of  high  ioniza¬ 
tion.  The  magnitude  of  the  polar  ion  imbalance  will  be  limited  by  the  ionic  space  charge  as 
expressed  by  the  next  statement. 

5.  The  effect  of  the  electric  field  created  by  an  imbalance  in  the  polar  ion  concentrations  limits  the 
ionic  imbalance  that  can  be  established  by  the  unequal  rates  of  diffusion  of  positive  and  negative 
ions  to  the  walls.  Since  the  electric  field  is  proportional  to  the  charge  enclosed  within  the 
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Fig.  1 1  —  Computed  and  measured  ion 

concentrations  at  distances  of  12  and  18  cm 
downwind  of  the  source.  The  upper  set  of  lines 
pertains  to  the  rectangular  channel;  the  lower  set 
of  lines  pertains  to  the  diffusion  cell  with  0.3-cm 
plate  spacing. 


volume  (Gauss’  Law),  the  limiting  effect  of  the  field  on  the  ion  ratio  is  strongly  dependent  on 
the  size  of  the  neutralizer. 

A  symmetric  charge  distribution  on  an  aerosol  does  not  imply  that  the  aerosol  is  in  charge 
equilibrium.  By  using  a  diffusion-dominated  region  downwind  of  an  ionizer,  it  is  possible  at 
some  point  to  balance  the  faster  diffusion  of  negative  ions  to  the  particles  by  a  predominance  of 
positive  ions  and  obtain  an  aerosol  with  zero  net  charge.  However,  this  is  a  transient  state,  not 
an  equilibrium  one,  and  should  not  necessarily  agree  with  the  predictions  of  an  equilibrium 
theory  of  charging. 

If  the  aerosol  concentration  is  much  greater  than  the  ion  concentration  (low  ionization  rate)  and 
if  the  net  charge  on  the  system  is  zero,  then  conservation  of  charge  requires  a  symmetric  charge 
on  the  aerosol.  This  case  is  usually  of  little  importance  in  aerosol  charging  because  the  time  for 
charging  at  low  ionization  rates  is  extremely  long.  However,  this  same  effect  can  play  a  signifi¬ 
cant  role  downwind  of  the  source  in  a  long  ionizer  with  larger  cross  section.  The  larger  cross 
section  prevents  a  large  ion  imbalance  (statement  5),  and  the  longer  residence  time  allows  the 
ion  concentrations  to  decay  to  levels  smaller  than  the  aerosol  concentrations.  Since  the  negative 
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ions  attach  more  readily  to  the  particles,  the  positive  to  negative  ion  ratio  increases  and  the  ratio 
of  charged  particles  approaches  unity  as  required  by  conservation  of  charge  (assuming  the  con¬ 
centrations  of  ions  are  negligible  compared  to  aerosols).  Examples  of  this  are  discussed  in  Sec¬ 
tion  5.1  and  in  Appendix  C. 


5.  EXPERIMENTAL  RESULTS 

5.1  Evaluation  of  NRL  Ionizers 

A  pair  of  our  early  ionizers  consist  of  sections  of  rectangular  waveguide  tubing  (aluminum)  with 
internal  dimensions  of  1.9  x  4.5  cm,  lengths  of  12.5  and  30  cm,  with  an  Am-241  foil  source  (256 
fiC i),  affixed  approximately  2  cm  from  the  inlet  of  each.  The  ends  are  sealed  by  plates  fitted  with 
tubing  adapters.  We  use  our  ionizers  in  conjunction  with  a  differential  mobility  analyzer  (DMA)  that 
has  a  sample  flow  rate  of  50  cc/s,  so  the  residence  time  of  aerosols  in  the  ionizers  is  roughly  2  and  5 
s,  respectively.  An  initially  perplexing  result  was  that  the  shorter  ionizer  consistently  produced  a  ratio 
of  negative  to  positive  particles  of  about  1.3,  whereas  the  longer  ionizer  produced  ratios  in  the  range 
of  1.8  to  2.1  (see  Fig.  12),  this  was  the  opposite  of  what  we  had  expected.  Our  initial  prediction  was 
that  in  the  longer  ionizer  the  ions  would  age  and  the  ratio  of  mobilities  would  decrease  as  observed  in 
laboratory  experiments,  thus  producing  a  more  balanced  charge.  However,  our  modeling  results  show 
that  the  shorter  ionizer  has  sufficient  numbers  of  ions  exiting  to  substantially  alter  the  aerosol  charge 
distribution  after  the  aerosol  has  left  the  ionizer.  The  small  diameter  of  the  tube  downstream  of  the 
ionizer  reduces  the  electric  field  due  to  the  ionic  imbalance  and  allows  the  negative  ions  to  diffuse  to 
walls,  relatively  unimpeded  by  the  predominance  of  positive  ions.  In  contrast,  the  longer  ionizer  has  a 
substantially  greater  residence  time  downwind  of  the  source,  allowing  the  ion  concentrations  to  decay 
to  a  much  lower  level  before  exiting.  The  electric  field  due  to  the  ionic  space  charge  prevents  a  large 
imbalance  in  the  polar  ion  ratio  while  in  the  ionizer  channel,  and  the  concentrations  of  the  exiting  ions 
are  now  too  low  to  substantially  modify  the  aerosol  charge  distribution  in  the  neutralizer. 

An  inverse  relation  between  aerosol  concentration  and  the  ratio  of  negative  to  positive  particles 
was  also  observed  when  the  30-cm-long  ionizer  was  used.  This  is  an  example  of  the  situation  dis¬ 
cussed  in  Section  2  and  Statement  7  of  Section  4,  where  the  particles  become  the  dominant  sink  for 
the  ions.  Figure  13  is  a  plot  of  the  charged  particle  ratio  vs  the  concentration  of  negative  particles 
within  a  single  channel  of  the  DMA.  For  the  specific  aerosol  used  during  this  experiment  (gas  phase 
nucleation  of  sulfuric  acid  and  ammonia  as  described  in  Appendix  A),  the  value  of  1000  particles/cc 
within  the  considered  interval  corresponds  to  a  size  distribution  (similar  to  that  presented  in  Fig.  A3) 
that  contains  on  the  order  of  20,000  particles/cc.  The  decrease  of  aerosol  charge  ratio  to  values  less 
than  1.5  required  total  aerosol  concentrations  greater  than  100,000/cc. 

5.2  Evaluation  of  Commercially  Available  Ionizers 

The  negative  to  positive  ratio  of  singly  charged  particles  exiting  two  ionizers  manufactured  by 
Thermo  Systems  Inc.  (TSI)  and  one  ionizer  produced  by  the  University  of  Minnesota  for  the  National 
Bureau  of  Standards  (NBS)  were  measured  as  a  function  of  particle  radius  by  using  the  aerosol  gen¬ 
eration  technique  described  in  Appendix  A  and  the  differential  mobility  analyzer  system  depicted  in 
Fig.  Al.  The  dominance  of  singly  charged  particles  was  ensured  by  the  steep  slope  of  the  size  distri¬ 
bution  as  shown  in  Figs.  A2  and  A3.  The  ionizers  tested  were  TSI  models  3012  and  3077,  and  the 
ionizer  loaned  to  us  by  NBS  (which  bore  no  identification).  The  behavior  of  the  TSI  model  3054  was 
also  investigated,  as  will  be  described  later.  All  of  these  ionizers  use  stainless  steel  encapsulated  Kr- 
85  (beta)  sources  ranging  from  2  to  10  millicuries  (mCi).  The  TSI  model  3012  is  fitted  with  screens 
interior  to  both  the  inlet  and  outlet,  a  feature  not  present  on  the  other  ionizers.  The  results  of  these 
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Fig.  12  —  Evaluation  of  the  particle  charge  ratio  produced 
by  two  early  NRL  ionizers  with  different  residence  times 


tests  are  presented  in  Fig.  14,  which  shows  that  none  of  these  ionizers  produced  a  truly  symmetric 
charge  distribution.  A  flow  rate  of  50  cc/s  was  used  throughout  the  experiment,  except  for  the  few 
points  noted.  The  physically  smallest  neutralizer,  TSI  model  3077,  provided  the  charge  distribution 
that  was  closest  to  symmetric  from  this  group  of  ionizers.  The  significantly  lower  particle  charge  ratio 
(N  _/N +)  obtained  with  the  TSI  3077  is  attributed  to  the  small-diameter  tubing  provided  for  inlet  and 
outlet  connections.  These  metal  tubes,  designed  to  reduce  radiation  leakage,  are  rather  long  and 
extend  into  the  active  region  of  the  ionizer,  so  that  the  tube  samples  high  ion  concentrations  that 
decay  during  passage  through  the  tube.  The  diffusional  loss  of  ions  within  the  small-diameter  tube 
results  in  an  abundance  of  positive  ions  that  alter  the  aerosol  charge  distribution  obtained  over  the 
source  (as  in  the  case  of  the  short  NRL  ionizer).  As  stated  earlier,  this  is  not  an  equilibrium  process 
and  the  results  will  depend  upon  the  flow  rate  and  aerosol  size  distribution.  In  an  effort  to  verify  this 
effect,  the  particle  charge  ratio  was  determined  for  a  TSI  model  3054  neutralizer  fitted  with  either  a 
7.0  or  3.5  mm  ID  metal  tube  32  cm  in  length,  at  the  exit.  Comparisons  were  made  between  cases  for 
which  the  metal  tubes  were  inserted  into  the  ionized  region  vs  those  for  which  the  tube  inlet  coincided 
with  the  normal  ionizer  outlet.  The  particle  charge  ratio  was  reduced  when  the  tube  was  inserted  into 
the  ionized  region  for  both  sizes  of  tube,  from  1.65  to  1.44  for  the  larger  diameter  tube,  and  from 
1.50  to  1.23  for  the  smaller  diameter  tube.  As  expected,  the  smaller  diameter  tube  permits  diffusion 
to  establish  a  larger  ion  ratio  before  it  is  limited  by  the  electric  field  (and  before  the  ion  concentra¬ 
tions  are  diminished  by  recombination).  The  fact  that  the  size  of  the  tube  connected  to  the  ionizer 
outlet  affected  the  particle  charge  ratio  indicates  that  sufficient  ions  exit  the  ionizer  to  niter  the  charge 
distribution  while  traversing  downwind  plumbing. 

We  emphasize  that  except  for  the  TSI  model  3077,  these  ionizers  are  designed  to  reduce  the  net 
charge  of  highly  charged  aerosols;  however  this  does  not  guarantee  zero  net  charge  or  charge 
equilibrium. 
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negatively  charged  particles  of  the  same  radius.  (The  total  concentration  of  particles  is  approximately  20  times  the  abscissa 
value.) 


6.  DESIGN  OF  AN  IONIZER  THAT  PRODUCES  A  NEAR-SYMMETRIC  EQUILIBRIUM 
CHARGE  DISTRIBUTION 

We  previously  stated  that  if  positive  and  negative  ions  have  different  mobilities,  there  can  be  no 
equilibrium  downwind  of  the  ion  source.  This  implies  that  removal  of  the  ions  immediately  downwind 
of  the  source  is  necessary  to  eliminate  the  wall  losses  that  alter  the  polar  ion  ratio.  A  neutralizer  must 
also  provide  sufficient  residence  time  for  the  aerosol  within  the  ionized  region  to  establish 
equilibrium.  Symmetry  of  the  charge  distribution  requires  that  the  ratio  of  positive  to  negative  ion 
concentrations  compensate  for  the  inequality  in  attachment  coefficients  of  positive  and  negative  ions. 
This  may  be  accomplished  by  incorporating  diffusion  surfaces  within  the  ionized  region  for  the  re¬ 
moval  of  more  negative  than  positive  ions;  these  surfaces  can  also  serve  to  reduce  the  electric  field 
associated  with  ion  imbalance.  The  ion  concentrations  must  be  low  enough  so  that  ion  diffusion 
losses  are  comparable  to  the  losses  due  to  recombination;  yet  the  ion  concentrations  must  also  be  sub¬ 
stantially  greater  than  the  aerosol  concentrations,  or  the  aerosol  charge  distribution  will  be  concentra¬ 
tion  dependent.  To  obtain  the  desired  trade-off  between  the  residence  time  required  for  aerosol 
charge  equilibrium  and  the  diffusive  loss  of  aerosols  to  the  walls,  it  is  necessary  to  tailor  the  ionizer 
to  a  specific  airflow  and  maximum  aerosol  load.  In  certain  instances  it  may  be  advantageous  to  bring 
the  aerosol  close  to  charge  equilibrium  with  a  front-end  ionizer  prior  to  the  final  charge  neutraliza¬ 
tion. 
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Fig.  14  —  The  ratio  of  negatively  to  positively  charged  particles 
obtained  with  four  different  ionizers 


An  ionizer  based  on  the  above  requirements  has  been  constructed  (see  Fig.  15),  using  a  Th-232 
sheet  4.0  cm  wide  and  18  cm  long  (total  activity  of  3.7  ^Ci)  in  conjunction  with  10  diffusion  plates 
that  are  perpendicular  to  the  active  sheet.  This  assembly  is  inserted  into  a  rectangular  section  of  tub¬ 
ing  (1.9  X  4.5  cm  cross  section  mentioned  earlier),  and  ion  removal  is  accomplished  by  a  4  x  3  cm 
plate  immediately  downwind  of  the  source  to  which  +20  V  is  applied.  As  shown  in  Fig.  14.  this 
NRL  prototype  ionizer  produced  the  most  symmetric  charge  distribution  of  the  ionizers  tested.  Esti¬ 
mates  of  the  polar  ion  concentrations  in  the  ionized  region  as  determined  from  measurements  immedi¬ 
ately  downwind  of  the  ionizer  indicated  positive  (negative)  ion  concentrations  of  40,000  (21,000)  per 
cc.  Consideration  of  the  aerosol  flow  rate,  ionizer  geometry,  and  ion  mobility  dictate  the  optimum 
voltage  and/or  length  of  the  ion  collection  plate  in  order  to  minimize  particle  loss.  The  minimum 
residence  time  within  the  ionizer  for  establishment  of  aerosol  charge  equilibrium  was  estimated  by 
using  the  aerosol  charging  model  discussed  in  Section  2.  The  diffusion  surface  area  was  determined 
by  trial  and  error  measurements  of  the  charged  particle  ratio.  The  loss  of  particles  larger  than  0.01 
/zm  in  this  ionizer  due  to  diffusion  is  not  severe;  both  measurements  and  calculations  indicate  that  the 
loss  of  0.01  nm  radius  particles  to  be  less  than  5%  at  the  prescribed  flow  rate  of  50  cm3  s_l. 
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Fig.  15  —  A  prototype  ionizer  that  produces 
a  near  symmetric  aerosol  charge  distribution 


7.  CONCLUSIONS 

The  charge  distributions  of  aerosols  exiting  radioactive  neutralizers  are  the  result  of  a  complex 
relationship  between  the  parameters  of  the  neutralizer  and  the  aerosol  charging  mechanism.  In  most 
(if  not  all)  of  the  present  ionizers,  the  charge  distribution  will  be  changing  downwind  of  the  source 
region  and  no  truly  equilibrium  distribution  can  be  maintained.  The  culprit  is  the  changing  polar  ion 
ratio  produced  by  the  combination  of  ionic  recombination  and  the  more  rapid  diffusion  of  negative 
ions  to  the  walls.  The  space  charge  electric  field  plays  an  important  role  in  limiting  the  degree  of 
ionic  imbalance  that  will  develop  through  unequal  diffusion  of  positive  and  negative  ions  to  the  walls. 
This  effect  is  strongly  influenced  by  the  size  of  the  ionizer. 

A  number  of  studies  (Hussin  et  al.  (1983),  Porstendorfer  et  al.  (1984),  Adachi  et  al.  (1985), 
Kousaka  et  al.  (1985))  have  been  initiated  to  test  various  charging  theories.  Each  of  these  studies  has 
relied  on  an  ionizer  to  bring  the  aerosol  to  “charge  equilibrium.”  This  study  suggests  that  these 
results  may  not  be  a  valid  test  of  the  theory,  since  the  measured  charge  distribution  is  most  likely  not 
only  determined  by  the  parameters  of  the  theory  but  also  by  the  particular  ionizer,  downwind  plumb¬ 
ing  geometry,  and  flow  rate. 

The  results  and  understanding  of  the  ionizer/neutralizer  provided  by  this  study  will  hopefully 
make  it  possible  to  develop  a  new  family  of  neutralizers  that  produce  a  well-defined  aerosol  charge 
distribution.  One  such  attempt  has  been  reported  in  this  report. 
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Appendix  A 

AMMONIUM  SULFATE  AEROSOL  GENERATION 
BY  GAS  PHASE  REACTION 


A  reproducible  and  stable  source  of  aerosol  particles  in  the  0.01  to  0.1  /on  size  range  has  been 
generated  by  the  following  technique.  Figure  A1  illustrates  the  aerosol  generation  and  measurement 
arrangement.  A  stream  of  charcoal  and  absolute  filtered  room  air  is  divided  into  three  streams:  acid, 
base,  and  dilution  flows.  The  acid  flow  (typically  1  cc/s  is  introduced  at  one  end  of  a  5-cm-ID  glass 
tube  containing  a  97%  H2S04  solution  that  has  a  surface  area  approximately  80  cm  long  and  4  cm 
wide.  The  acid  is  contained  behind  a  glass  dam  at  the  downwind  end  of  the  surface.  Air  flowing 
over  the  acid  contains  very  low  concentrations  of  sulfuric  acid  vapor.  The  base  flow  (also  typically  1 
cc/s)  is  passed  over  the  surface  of  a  0.05%  solution  NH4OH  contained  in  a  1-liter  flask  and  exits  from 
the  top.  The  air  containing  NH3  is  combined  with  the  air  containing  acid  vapor  immediately 
downwind  of  the  dam.  Aerosol  formation  by  binary  homogeneous  nucleation  is  almost  immediate. 
The  flows  are  allowed  to  mix  for  the  remaining  length  of  the  glass  tube  (40  cm)  and  then  combined 
with  the  dilution  flow  (typically  48  cc/s)  and  used  for  the  experiment.  The  use  of  a  large  volume 
holding  chamber  is  not  recommended  (unless  adequate  mixing  is  available),  for  it  was  found  that  the 
aerosol  concentrations  were  much  more  stable  without  it.  Examples  of  aerosol  size  distributions  are 
presented  in  Figs.  A2  and  A3  for  which  the  airflows  over  both  he  acid  and  the  base  solutions  were 
increased  from  1  to  5  cc/s.  As  shown  in  these  figures  the  aerosol  peak  is  shifted  toward  smaller  sizes 
when  the  residence  times  of  the  airflows  over  their  respective  solutions  is  decreased.  The  reduction 
of  the  surface  area  of  the  acid  bath  also  decreases  the  mean  particle  size.  An  aerosol  of  larger  size 
particles  can  be  produced  by  the  redirecting  the  airflow  that  has  passed  over  the  base  solution,  over 
the  acid,  producing  size  distributions  as  shown  in  Fig.  A4. 
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Fig.  A I  —  Aerosol  generation  and  measurement  apparatus 
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Typical  size  distribution  for  1  cc/s  flows 
over  acid  and  base  solutions 
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Fig.  A3  —  Typical  size  distribution  for  5  cc/s  flows 
over  acid  and  base  solutions 


Appendix  B 

MEASUREMENTS  OF  ION  CURRENT  VS  COLLECTOR  VOLTAGE 


Figures  B1  through  B4  are  I-V  curves  obtained  from  the  experimental  arrangement  depicted  in 
Fig.  10,  where  the  radioactive  sources  consisted  of  a  1.0  cm  diameter,  5.3  /*Ci  Am-241  foil.  Five 
grounded  vertical  shield  wires  were  placed  just  upwind  of  the  plate  in  an  effort  to  reduce  the  upw  ind 
penetration  of  the  field  from  the  plate.  Figure  B1  shows  that  when  the  distance  between  the  foil  and 
collector  was  2  cm,  the  linear  (conduction)  region  of  the  I-V  curve  extends  from  about  -2  to  +2 
volts,  and  the  curve  flattens  as  the  majority  of  the  ions  are  collected.  The  slope  of  the  curve  does  not 
go  to  zero  because  increased  voltage  reduces  the  ion  loss  by  recombination;  this  also  applies  to  Figs. 
B2  and  B3  for  which  the  source  to  collector  distances  were  6.0  and  30  cm  respectively.  The  further 
the  ion  collector  is  from  the  source,  the  easier  it  is  obtain  current  saturation  above  2  V  because  the 
rate  of  recombination  in  the  ion  collector  region  is  reduced.  The  numbers  near  the  I-V  curves  give 
the  slope  of  the  linear  portion  of  the  curve  in  units  of  Amp  Volt’  '. 

The  relative  mobilities  of  the  ions  can  be  estimated  from  the  the  slopes  of  the  conductive  region 
and  the  total  ion  concentration  by  the  use  of  the  following  relation: 


k_  X_n, 


where 

k  are  the  ion  mobilities, 

X  are  the  slopes  (polar  electrical  conductivities),  and 
n  are  the  ion  concentration. 

Evaluation  of  this  relation  from  Fig.  B3  (d  =  30  cm)  resulted  in  a  ratio  of  negative  to  positive 
ion  mobilities  of  approximately  1.4.  This  is  significantly  greater  than  the  ratio  of  ion  mobilities  we 
used  in  our  aerosol  charging  theory  (Hoppel  and  Frick.  1986).  which  are  applicable  for  atmospheric 
ions.  The  same  theory  was  used  in  Section  2  of  this  report:  however  it  was  revised  to  accommodate 
the  greater  difference  in  ion  mobilities.  Our  calculations  of  ion  deposition  in  a  channel,  given  in  Sec¬ 
tion  4,  used  mobilities  that  gave  a  ratio  of  1 .4. 

Figure  B4  illustrates  the  large  ratio  of  positive  to  negative  ions  that  can  be  obtained  with  a  diffu¬ 
sion  cell.  For  this  case  the  radioactive  source  was  located  12  cm  from  the  shield  wires,  and  the  6  cm 
immediately  in  front  of  the  shield  wires  was  occupied  by  a  "diffusion  cell"  (a  set  of  five  parallel 
plates,  each  4.0  x  6.0  X  0.05  cm,  discussed  in  Section  3.2). 
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Appendix  C 

TIME-DEPENDENT  CHARGING  CALCULATIONS 


We  have  modeled  the  time-dependent  charging  of  a  monodisperse  aerosol  by  applying  a  Runge- 
Kutta  method  to  the  combined  set  of  Eqs.  (5)  to  (11).  Particles  of  radius  0.01  nm  were  used  in  all 
but  one  of  the  following  figures  because  the  0.01-/tm  particle  is  within  the  size  range  of  our  interests 
and  lacks  the  multiple  charges  that  tend  to  clutter  the  figure.  The  concentration  of  10,000  cm-3  was 
chosen  as  a  standard  aerosol  concentration  for  comparisons  between  the  figures,  thus  the  results 
presented  here  would  not  apply  to  concentrations  much  greater.  The  bipolar  ion  production  rate  of 
100,000  ion  pairs  cm-3  was  selected  as  representative  of  an  instance  for  which  the  recombination  loss 
of  ions  may  not  necessarily  be  the  predominant  loss  mechanism.  Although  wall  losses  are  not 
included  in  the  aerosol  charging  calculations,  we  are  considering  regimes  for  which  ion  losses  to  walls 
would  appreciably  alter  the  balance  of  ions.  (The  effect  of  walls  on  the  ratio  of  the  ion  concentra¬ 
tions  was  simulated  by  using  different  initial  ion  concentrations,  as  shown  in  Fig.  8.) 

In  Figs.  Cl  to  C4,  concentrations  of  positive  (negative)  ions  are  identified  by  n  +  («_),  the 
uncharged  particles  by  Os,  and  the  charged  particles  by  the  appropriate  sign  followed  by  the  number 
of  charges. 

Figure  Cl  shows  the  charging  of  0.01-jnm  particles  for  an  ionization  rate  of  106  ion  pairs/cm3/s. 
The  relaxation  time  is  reduced  to  ~3  s  compared  to  10  s  for  Fig.  6  (105  ion  pairs/cm3/s). 

Figure  C2  shows  that  if  there  were  preexisting  5  X  105  concentrations  of  bipolar  ions,  the 
relaxation  time  is  only  reduced  to  about  9  s. 

In  Fig.  C3  the  entire  population  of  105  particles  has  been  given  a  negative  single  charge,  and  the 
relaxation  time  is  about  15  s;  note  the  equilibrium  values  for  positive  and  negative  particles.  Figure 
C4  shows  that  if  an  initial  ion  imbalance  is  introduced  so  as  to  eliminate  the  net  space  charge,  the 
charge  equilibrium  values  will  be  the  same  as  those  in  previous  figures.  This  difference  in 
equilibrium  is  due  to  the  charge  conservation  built  into  the  governing  equations  and  results  in  the 
majority  of  the  space  charge  being  accommodated  by  an  imbalance  in  the  equilibrium  ion  concentra¬ 
tions,  with  the  remainder  carried  on  the  particles. 
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Appendix  D 

ION  AEROSOL  ATTACHMENT  COEFFICIENT  CALCULATION 


It  is  clear  from  the  calculations  carried  out  in  this  report  that  the  physical  properties  of  the  ions 
are  important  factors  in  determining  the  ion-aerosol  attachment  coefficients.  Particularly  important  is 
the  ratio  of  the  positive  to  negative  diffusion  coefficient  in  establishing  the  ratio  of  positively  to  nega¬ 
tively  charged  particles  (Eqs.  (9)  to  (11)).  The  ion  species  changes  with  age  as  a  result  of  charge 
exchange  between  molecules  of  different  ionization  potential  or  electron  affinities,  and  chemical  and 
clustering  processes.  Ionic  mass  as  a  function  of  age  has  been  measured  by  a  number  of  investiga¬ 
tors,  but  the  ionic  chemistry  is  extremely  complex  and  depends  on  the  concentration  and  types  of 
trace  gases  in  the  air  [see  for  example,  Huertas  and  Fontan  (1975);  Huertas  et  al.  (1978);  Mohnen 
(1977);  and  Perkins  and  Eisele  (1984)].  There  is  currently  no  consensus  on  the  ionic  species  that 
exists  in  “real”  air  at  atmospheric  pressure  on  the  time  scales  of  interest  in  the  neutralizer.  Most 
definitive  data  have  come  from  mass  spectroscopic  measurements  and  at  times  that  correspond  to 
microseconds  when  converted  to  atmospheric  pressure.  Since  the  molecular  collision  rate  is  of  the 
order  of  1010  s-1,  the  ion  will  suffer  collisions  with  trace  gases  at  the  ppb  level  in  less  than  a  second. 
The  chemistry  of  the  initial  ion  probably  bears  little  resemblance  to  the  ions  that  exist  after  tens  of 
milliseconds.  Measured  values  of  ion  mobilities  (and  recombination  coefficients)  seem  to  fall  into 
tnree  classes  of  measurements;  (1)  initially  formed  ions  of  very  short  lifetime  made  in  ultrapure  air  at 
low  pressures,  (2)  ions  with  a  lifetime  of  the  order  of  a  second  in  artificially  ionized  air  at  about 
atmospheric  pressure,  and  (3)  atmospheric  ions  that  have  lifetimes  of  hundreds  of  seconds  and  exist  in 
concentrations  of  100  to  1000  cm"3.  Atmospheric  ions  have  reduced  mobilities  of  about  1.15  and 
1.25  cm2  V-1  s'1  for  positive  and  negative  ions  respectively  [Mohnen,  (1977)]  and  a  recombination 
coefficient  of  1.4  x  I0~6cm3s_1  [Nolan,  (1943)]. 

There  have  been  many  measurements  of  ion  mobilities  in  the  laboratory,  with  ions  generated  by 
various  artificial  ionization  sources  dating  back  to  the  early  1900s  when  ions  were  first  discovered. 
In  these  experiments  the  ion  concentrations  are  much  greater  than  in  the  atmosphere,  and  therefore 
the  ion  lifetime  is  much  less.  For  instance,  the  recombination  lifetime  of  ions  at  a  concentration  of 
106  is  less  than  1  s.  A  table  of  mobility  values  measured  by  different  authors  in  air  of  various 
degrees  of  purity  and  various  ionization  sources  is  given  by  Mohnen  (1977).  The  average  mobilities 
of  these  ten  measurements  are  1.33  and  1.84  cm2  V-1  s-1  at  STP  with  a  negative  to  positive  mobility 
ratio  of  about  1.4.  The  conditions  under  which  these  measurements  were  made  are  more  like  those 
that  must  exist  in  ionizers  than  atmospheric  conditions. 

In  our  earlier  work  [Hoppel,  (1985);  Hoppel  and  Frick  (1986)];  we  were  interested  primarily  in 
atmospheric  ions  and  we  used  positive  and  negative  mobility  values  of  1.2  and  1.35  cm2  V-1  s~'. 
This  gives  a  negative  to  positive  mobility  ratio  of  only  1.12.  While  these  values  are  probably  the  best 
available  for  aged  tropospheric  ions,  they  are  clearly  not  the  best  values  for  conditions  which  exist  in 
ionizers.  We  have  therefore  recalculated  all  the  ion-aerosol  coefficients  using  positive  and  negative 
mobility  values  of  1.3  and  1.8,  respectively. 
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In  addition  to  the  mobilities,  the  ionic  masses  are  also  required  for  the  calculation  of  the  ionic 
mean-free-paths.  The  mean  free  path  is  directly  proportional  to  the  mobility  and  inversely  propor¬ 
tional  to  the  factor 

r  ''i  1/2 


which  is  a  much  weaker  dependence  since  the  ionic  mass  Mion  is  much  greater  than  the  mass  of  the 
air  molecule  Mair  (28.8  AMU).  There  is  not  a  unique  relationship  between  mobility  and  ionic  mass; 
but  empirical  plots  of  mass  vs  mobility  [Cabane  et  al.  (1976),  Huertas  et  al.  (1974),  and  Meyerott  et 
al.  (1980)]  indicate  that  for  the  mobilities  used  here,  ionic  masses  of  200  AMU  and  100  AMU  for 
positive  and  negative  ions  are  very  reasonable  values. 

The  ion-aerosol  attachment  coefficients  have  been  recalculated  for  positive  and  negative  ions 
with  the  properties  cited  above  and  are  given  in  Tables  D1  and  D2  and  plotted  in  Fig.  Dl.  The  pro¬ 
cedure  used  is  that  presented  in  Hoppel  and  Frick  (1986)  with  one  modification.  In  the  case  of  non- 
symmetric  ion  properties,  the  electrostatic  energy  gained  by  the  two  ions  during  an  encounter  is  dif¬ 
ferent,  and  the  energy  loss  required  for  trapping  is  different  for  the  two  ions.  Instead  of  adjusting 
the  ion-aerosol  trapping  distance  using  the  functional  relationship  given  in  our  earlier  paper,  we  here 
adjusted  the  energy  required  to  ensure  trapping.  The  energy  loss  required  to  ensure  trapping,  calcu¬ 
lated  for  ion-ion  recombination  was  then  partitioned  between  the  two  ions  (for  the  ion-aerosol 
encounter)  in  a  manner  which  accounts  for  the  fact  that  the  negative  (lighter)  ion  acquires  more  elec¬ 
trostatic  energy  than  the  positive  ion  and  the  fact  that  the  negative  ion  loses  more  energy  during  a  col¬ 
lision  with  a  neutral  molecule  because  it  is  lighter.  The  ion-aerosol  trapping  distances  are  then  calcu¬ 
lated  from  Eq.  (34)  of  Hoppel  and  Frick  (1986)  separately  by  using  different  energies  for  positive  and 
negative  ions.  Both  the  old  and  new  method  give  the  same  result  when  the  particles  are  small  (r  < 
0.02  /;m),  but  when  the  polar  ions  have  large  differences  in  mobilities  as  is  the  case  here,  then  the  old 
method  leads  to  errors  at  intermediate  sizes. 
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Table  D1  —  Attachment  Coefficients  For  200  AMU  Ions  (cm3  s  ')  (multiply  by  106) 
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